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Using HF+BGS method we study light nuclei with nuclear charge in the range 2 < Z < 
8 and lying near the neutron drip line. The HE method uses effective Skyrme forces 
and allows for axial deformations. We find that the neutron drip line forms stability 
peninsulas at s^He and These isotopes are found to be stable against one neutron 
emission and possess the highest known neutron to proton ratio in stable nuclei. 


1. Introduction 

Nuclei with extreme neutron access lying close to the neutron drip line constantly 
attract interest both in theoretical and experimental research. Experiments with 
radioactive nuclear beams that are conducted in Dubna, GSI, RIKEN have opened 
new opportunities for obtaining exotic nuclei with impressive neutron excess. The 
construction of FRIB in MichigaiP (the operational run is planned in 2022) would 
boost experimental abilities to approach the neutron drip line. Among major fun¬ 
damental microscopic approaches in the studies of nuclei with neutron excess one 
can name the Hartree-Fock-Bogoliubov approach (HFB), Hartree-Fock plus BCS 
pairing (HF-fBCS) with effective forces and relativistic mean field theory.l^lt^ 

In Refs. [HI m Uni HU o im [m [m we explored the possibility of existence of 
islands and peninsulas of stability of nuclei with large neutron excess lying beyond 
the conventional neutron drip line. The calculations in these papers were done with 
the HF method using various type of effective Skyrme forceJ^^^^ and accounting 
for the axial deformation; the pairing was treated in the BCS scheme. We have 
demonstrated that in the regions of the nuclear chart corresponding to extreme 
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neutrons excess around magic and “new” magic numbers N = 32, 58, 82, 126, 184, 
258 there may exist peninsulas of neutron stability stretching beyond the conven¬ 
tional neutron drip line. In this picture a neutron rich nucleus, which is unstable 
against neutron separation, may regain stability if one adds to it a certain num¬ 
ber of neutrons thereby shifting it to the stability peninsula. In the HF approach 
this stability restoration is a result of complete filling of neutron subshells with 
large angular momentum. These subshells have a large centrifugal barrier and be¬ 
ing partially filled they are located in the continuous spectrum, which makes the 
corresponding nuclei unstable against neutron emission. However, when the neu¬ 
tron number increases these subshells descend into discrete spectrum and nuclei 
regain stability. The main aim of the present paper is to find stable nuclei with the 
highest possible neutron to proton ratio N / Z, which we believe sets a theoretical 
upper bound for this ratio. We shall consider only nuclei with 2 < Z < 8, since 
it has been shown earlierthe maximal N / Z ratio can be attained only for 
light nuclei. The largest values of N are obtained at stability peninsulas. The most 
extended stability peninsulas are obtained with the set of Skyrme forces, 

therefore, we shall use mostly the Ski Skyrme parameters. Let us note that various 
types of Ski forces were obtained and effectively used in Ref. [20l The same type of 
forces was used in Ref. [22] for the description of light exotic nuclei, however, without 
accounting for deformation, which is, certainly, a disadvantage compromising the 
predictive power of the obtained results. 


2. Methods and Results 

The detailed exposition of the method of solving the deformed Hartree-Fock (DHF) 
equations is given in Refs. |9| |23l HU |25l Here we present the results of the calcu¬ 
lations using SkI2 set of Skyrme forces. One particle wave functions in DHF are 
expanded in the axially deformed harmonic oscillator basis. The principle quantum 
number for the harmonic oscillator basis usually does not exceed A^o = 18 (amount¬ 
ing to 1330 basis functions). This dimension of the basis is more than necessary for 
the calculations of the isotopes with 2 < Z < 8, which provides high accuracy and 
reliability of obtained results. 

In case of nuclei lying on stability peninsulas along with the DHF calculations 
we do additional calculations, where HF equations are solved directly in coor¬ 
dinate space under assumption of spherical symmetry (SHF method) This is 
reasonable because nuclei belonging to stability peninsulas indeed possess spher¬ 
ical symmetry.l^E^ SHF calculations make possible the analysis of the role of 
continuum states in the cases when the DHF solution is spherically symmetric. 
We used the BCS constant pairing, where the pairing constant is set equal to 
G = {19.5/{N + Z))[l ± 0.5I(A^ - Z)/{N + Z)]^ “-b” and correspond to pro¬ 
tons and neutrons respectively. In DHF calculations pairing is restricted to bound 
one-particle states. This choice of pairing as well as the role of continuum states in 
the structure of nuclei lying at stability peninsulas is discussed in Refs. [HI [121 [13 
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[T1[T51 Let us note that the inclusion of continuum states into present SHF calcu¬ 
lations of nuclei with 2 < Z < 8, which form stability peninsula, did not affect the 
results that were obtained without inclusion of continuum states. This is, however, 
typical for magic and quasi-magic nuclei. Similar behavior has been observed for 
stability peninsulas in other parts of the nuclear chart 

Fig. 1 shows nuclear chart, where the squares represent nuclei that are sta¬ 
ble against one neutron emission in DHF calculations with Skl2 forces. Exactly the 
same picture is obtained with Skll forces. Grey squares on this chart are experimen¬ 
tally known stable nuclei. We determine one neutron drip line from the condition 
S'ra = 0, where S'nO denotes one neutron separation energy. Thereby, any positive 
value of Sn means stability against one neutron emission, even if S'„ is marginally 
small. We calculate one neutron separation energies assuming the validity of ap¬ 
proximation made in Koopmans theorem. Eor given Z the position of the neutron 
drip line is determined as follows. We start with a stable isotope and increase the 
neutron number until stability is lost and one neutron separation energy changes 
its sign. Knowing that one-neutron drip line can exhibit nonmonotonic behavior^ 
we continue increasing N. If at some larger value of N the value of Sn changes 
its sign again then this means that we have located a stability peninsula. As we 
have mentioned such effect of stability restoration appears when neutron subshells 
with large angular momentum getting fully filled intrude from continuous spec¬ 
trum into the bound spectrum.EUi^ (IiP^ the formation of stability peninsulas was 
alternatively termed irregular behavior of the neutron drip line). Using SkI2 set 
of Skyrme forces we have found one stability peninsula situated at ^®He with the 
fully filled subshell ld 5/2 and another one at with the fully filled subshell lf 7 / 2 . 
These isotopes that disrupt the monotonic behavior of the drip line are illustrated 
in Fig. 1, where they are highlighted with red (color online) color and have arrows 
pointing at them. Both of these nuclei possess spherical symmetry, which is typical 
for nuclei with fully filled shells. The neutron to proton ratio in ^®He and is 
N/Z = 8 and N/Z ~ 5.67 respectively, which is substantially larger than the so far 
known neutron to proton ratios for stable nuclei.l^®^ Fig.l also shows the stability 
peninsula consisting of one isotope Later we shall discuss in detail the stability 
restoration for this isotope. 

Fig. 2 shows one neutron separation energies as functions of A = {Z -\- N) ob¬ 
tained with SkI2 Skyrme forces. One can clearly see S'„ changing sign at ^®He, which 
results from completely filled subshell ld 5/2 being bound. One neutron separation 
energy of this isotope as estimated according to the Koopmans theorem equals 0.45 
MeV. Let us consider HF potentials corresponding to the state ld 5 / 2 , where this 
state is weakly bound. Because the nucleus is spherical we can use SHF method, 
which has the advantage of representing potentials and wave functions directly in 
coordinate space.^^ The DHF value Sn — 0.450 is close to the SHF value Sn = 
0.406 MeV. Fig. 3 shows the Fermi level ld 5/2 for ^®He and the corresponding SHF 
potential. The dotted line depicts the wave function of the Fermi level, which is 
localized under the centrifugal barrier. The height of the centrifugal barrier in the 
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HE potential equals 1.35 MeV, which should additionally enhance the stability of 
^®He against one neutrons emission. When we apply the BCS scheme in case of SHF 
calculations we also include those localized quasibound states with positive energy, 
whose wave functions are localized under the centrifugal barrier. Other continuum 
states are not considered. In case of ^®He and (Skl2 forces) we considered all 
localized quasibound states. The states were discretized by introducing boundary 
conditions, which make wave functions vanish outside the sphere with the radius 
40 fm. We found that the pairing remained zero even when we included localized 
quasibound states. 

Fig. 4 shows one neutron separation energies for carbon calculated with Skl2 
forces using the DHF method. One can see that Sn increases for A = 26, which 
signals the magicity of the number TV = 20 near the drip line. The neutron subshell 
lf 7/2 submerses into bound spectrum becoming completely filled for N = 34, which 
provides stability against one neutron emission for At this point Sn = 0.597 

MeV in DHF calculations and Sn = 0.594 MeV in SHF calculations. is a 
spherically symmetric nucleus. Similar to the case of ^®He in Fig. 5 we consider the 
SHF potential corresponding to the Fermi level lf 7/2 for The SHF potential 
for this state has a centrifugal barrier with the height 2.48 MeV, which additionally 
enhances the stability of this isotope. In Fig. 4 one witnesses the stability restoration 
for (26 neutrons) , which is again connected with the state lf 7 / 2 . Unlike the 
nucleus has quadrupole deformation /3 = —0.32, and its last filled neutron level 
has quantum numbers D = 7/2“. This level, which is the member of the multiplet 
lf 7 / 2 ) intrudes into the bound part of the spectrum when N = 26. 

Let us consider in more detail the isotope chains of He and C, which pass 
through ^®He and In Figs. 6, 7 one can see the neutron and proton root 
mean square (rms) radii for helium and carbon. The arrows indicate the 

values for ^®He and which form stability peninsulas at V = 16 and N = 34. 
Generally in neutron rich nuclei the difference of neutron and proton rms radii 
^Rn,p = may take large values and one can speak about neutron 

skin or neutron halo effects when Ai?„,p > 0.2 fm.I^ From the figure one can see 
large values of Ai?„_p for He and C isotopes in particular for ^®He and '^°C, which 
is a manifestation of a large neutron halo in these nuclei. In DHF calculations for 
^®He we obtained Ai?„^p = 4.288 — 2.150 = 2.138 fm and ~ 1.994. 

For 40c we got Ai?„,p = 4.709 - 2.800 = 1.909 fm and )U2/^p2^i/2 ^ ^ 532 . It 
is useful to compare these values with those of extremely neutron rich ^o 0,1^ where 
Ai?„,p = 1.29 fm and — 1.436 (SkM* forces). This suggests that 

there is a very large neutron halo in 4®He and ^Oc in the picture corresponding to 
SkI2 forces. 

Figs. 8, 9 show neutron and proton density distributions of 4®He and ^Oc calcu¬ 
lated within SHF approximation using SkI2 forces. SHF approximation yields the 
wave function in coordinate representation and is thus more appropriate for graph¬ 
ical purposes. Its use in this case is justified by the spherical form of these nuclei. 
Neutron halo in neutron rich nuclei is first of all characterized by the spatially ex- 
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tended neutron density distribution dominating over the “normal” proton density 
distribution. In this sense Eigs. 8, 9 provide a more clear picture of neutron halo 
in nuclei forming stability peninsulas. Neutron and proton density distributions in 
Figs. 8, 9 are typical for neutron rich nuclei.^ Our analysis shows that the tail of 
neutron distribution results from weakly bound one particle states. 

3. Conclusion 

We did HF + BCS calculations of neutron rich nuclei with charges 2 < Z < 8. 
We have demonstrated that beyond conventional neutron drip line there may exist 
stability peninsulas formed by ^®He and These nuclei set the record for the 
largest so far known neutron to proton ratio for stable nuclei. They are spherical 
in shape and have a spatially extended neutron halo. These newly found stability 
peninsulas complement our previous findings and their existence agrees with the 
general nature of the neutron drip line forming peninsulas at magic and quasimagic 
neutron numbers. 
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Fig. 1. (Color online). Fragment of NZ chart of even-even nuclei. Empty squares are nuclei sta¬ 
ble against one neutron emission in DHF calculations with SkI2 forces. Grey squares represent 
experimentally known stable nuclei. Black squares are stability peninsulas. 



Fig. 2. One neutron separation energies Sn of helium isotopes as a function of the neutron number 
obtained in DHF calculations with SkI2 forces. The arrow points at ^®He. 
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Fig. 3. The last filled level ld 5/2 ^nd the corresponding SHF potential V{r) for ^^He (SkI2 forces). 
The energy is 0.406 MeV, the centrifugal barrier has the height 1.35 MeV. The dotted line is the 
wave function of the level Idc/o. 
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Fig. 5. The last filled level lf 7/2 Ihe corresponding SHF potential V{r) for (SkI2 forces). 
The energy is 0.594 MeV, the centrifugal barrier has the height 2.48 MeV. The dotted line is the 
wave function of the level lf 7 / 2 - 
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Fig. 6. Neutron (filled squares) and proton (empty squares) root mean square radii of helium 
isotopes as functions of mass number. These are results of DHF calculations with SkI2 forces. 
Half-filled squares correspond to nuclei that are unstable against one neutron emission. 
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Fig. 8. Neutron and proton density distributions for ^®He calculated within SHF approach with 
SkI2 forces. 
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Fig. 9. The same as Fig. 8 but for 




